MG13 Proceedings: A lattice Universe as a 
toy-model for inhomogeneous cosmology 



! J- -P- Bruneton 1 

^amur Center for Complex systems (naXys), University of Namur, Belgium 
E-mail: jpbr@math.fundp.ac.be 



o 

(N 



Abstract. We briefly report on a previously found new, approximate, solution to 
Einstein field equations, describing a cubic lattice of spherical masses. This model 
mimics in a satisfactory way a Universe which can be strongly inhomogeneous at small 
scales, but quite homogeneous at large ones. As a consequence of field equations, the 
lattice Universe is found to expand or contract in the same way as the solution of a 
Friedmann Universe filled with dust having the same average density. The study of 



observables indicates however the possible existence of a fitting problem, i.e. the fact 
that the Friedmann model obtained from past-lightcone observables does not match 
with the one obtained by smoothing the matter content of the Universe. 
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CO ! 1. A toy-model for inhomogeneous cosmology 
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As a toy- model of inhomogeneous Universe, we study a lattice Universe of, e.g., galaxy- 
like objects. The parameters of the lattice are the comoving size L of the cells, the mass 
M inside each cell, and 77 the (comoving) "size" of the masses themselves. This means 
that we choose to describe the source term by a three dimensional comb of Gaussians 
of width 77. Two dimensionless ratios then govern the physics of the model: the typical 
depth of gravitational potential, GM/Lc 2 , and the degree of inhomogeneity of the lattice 
Universe, rj/L. When r]/L is of order 1, the lattice Universe is nearly homogeneous at 
all scales with very low density contrasts while small values rj/L <^ 1 corresponds to a 
Universe essentially made of voids and of small regions of very dense matter. Relevant 
parameters, phenomenologically speaking, are typically galaxy-sized masses, M ~ 10 11 
Solar masses, and L ~ IMpc, for which the ratio GM/Lc 2 is very small, of order 1CU 8 . 

This last point suggests to try and solve the Einstein equations perturbatively, as a 
power series in GM/Lc 2 . Doing so, Einstein equations are effectively linearized up to 
order three in the parameter y/M/L which turns out to be the relevant quantity for 
Taylor expanding the metric of spacetime and the stress-energy tensor Moreover, 
one can take advantage of the periodicity of the lattice, and thus look for a periodic 
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solution. For simplicity, we have assumed flat space with no cosmo logical constant. The 
solution then reads [I], in synchroneous gauge, g 00 — ~ 1> 9oi — 0, and 
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where £ = ±1 corresponds to an expanding or contracting lattice, and the function 
which accounts for the anisotropies of the gravitational field reads: 
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This metric solves Einstein equations G M „ = 87rG / 'c with no cosmological constant, 
and for a source term representing the lattice up to order (M/L) 3//2 , namely a stress 
energy tensor whose only non-vanishing component reads 
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where = t)^ + y^-hfy + 0(M/L) defines h$ via Eq. ([T]). The solution is typically 

valid for a limited range in time and space At and Ax, with Aa; ~ At <C Ly^L/M ~ 
lOGpc. 



(i) 



The metric enables one to compute the effective scale factor for such a Universe, and 
thus an effective Hubble parameter. Up to the order of approximation considered, one 
finds pp 



H(t) = H - *-Hlt + ^Hlt 



0{H%), (4) 
is identified with the comoving average density 



where the Hubble parameter at time t 
Hq = e \J~^~\J ^jjr- This corresponds to a flat Friedmann-Lemaitre model with equation 
of state w = 0. Thus, the model with discrete masses on a cubic lattice is identical to 
a Friedmann model with dust, with the corresponding energy density. This means that 
one cannot distinguish the distribution of mass (localized or smeared in an homogeneous 
manner) from purely kinematical considerations. In the language introduced in [2] , the 
model exhibits no strong backreaction, ie. the scale factor of the averaged Universe is 
identical to the average of the scale factor, up to the order of approximation considered, 
in accordance with other studies in similar contexts, see [3j H] . 
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2. Observables and the fitting problem 

This sort of commutativity property is however broken at the level of observables. In 
[5], we have studied the propagation of bundle of light rays in the spacetime given by 
the previous solution. By solving the geodesic equation, we have computed, up to order 
GM/Lc 2 , the redshift as a function of the affine parameter A, while solving perturba- 
tively the Sachs equations governing the expansion and the shear of the bundle, we have 
derived an analytical expression for the angular distance as a function of A. Combining 
both, this leads to distances as a function of redshift. 

Doing so, one finds that observables are given by their standard expression in the analog 
Friedmann model with the average density, plus non-trivial corrections that need not 
to stay negligible in general^]. This indicates that a weak backreaction [2], or a fitting 
problem [BJ, is present in the model. By carefully analysing the analytical expressions 
found, we have been able to establish a quantitative criterion, connecting the three 
length scales of the problem, that decides whether deviations to the laws of homoge- 
neous cosmology are to stay small. More precisely, we showed that in a Universe where 
masses are smeared enough, namely, satisfying 



the perturbations to Friedmann's expression for the observables remain small. A viola- 
tion of this bound, on the other hand, ie. when objects are too compact, or the Universe 
too inhomogeneous, implies a breakdown of the perturbative expansion, as far as the 
calculation of observables are concerned, although the pertubative scheme remains rel- 
evant in order to estimate the spacetime geometry. The next step is therefore to study 
more precisely the behavior of observables by solving non-perturbatively the coupled 
system of Sachs equations in the perturbed metric Eq. (CD). Indeed, the a/ M/L pertur- 
bative solution to Sachs equations cannot be trusted anymore in this case, because it 
decouples the shear from the expansion of the bundle of light rays, whereas the shear 
is expected to be large, and should thus play an important role. This is work in progress. 

In conclusion, we have found a realistic model for inhomogeneous cosmology that ex- 
hibits no strong backreaction, but weak backreaction whose amplitude is linked to the 
Weyl tensor sourcing for the shear, and therefore also linked to the "compactness" of 
the sources, or, if one wishes, to the degree of inhomogeneity of the Universe. 
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